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Abstract 

Explicit expressions for anharmonic contributions to the thermodynamic prop- 
erties with allowance for higher-order phonon-phonon interactions for closed-packed 
crystals are given, and the calculations for some fee metals near the melting (Ir, 
Rh) and martensite phase transition (Ca, Sr) points are carried out. A detailed 
comparison of anharmonic and electron contributions to the heat capacity of these 
metals is carried out. The computational results for high-temperature heat capacity 
agree well with the available experimental data. 
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Anharmonic effects (AE) in the lattice dynamics and thermodynamics have been ex- 
tensively discussed in the literature [1-5]. However the quantitative information on AE, 
both experimental and theoretical, is rather poor so far. The experimental difficulties are 
due to the fact that at high temperatures (T ~ T m , where T m is the melting temperature) 
when the AE in the thermodynamic properties become noticeable, it is usually difficult 
to identify their contribution to the heat capacity and thermal expansion from the con- 
tributions of vacancies and other thermally excited lattice defects [2,5-8]. The theoretical 
calculations, on the other hand, are based, as a rule, on rough models and, therefore, do 
not give any definite information on the AE in thermodynamic of real crystals. At the 
same time this information is needed for understanding the high temperature properties of 
solids. For example, it is important to clear up what is the role of self-anharmonic effects 
associated with phonon-phonon interactions, and what is the role of quasiharmonic ones 
which are only due to the change in the volume because of thermal expansion. Simple 
estimations M show that the AE contribution to heat capacity at high temperatures of 
the order of the electron contribution to heat capacity and has the same (linear) tem- 
perature dependence; the real ratio of these two contributions essentially depend on the 
features of the electron structure near the Fermi level Ep. However this problem has 
not been studied yet. All these problems are of special importance for refractory metals 
(T m > 2000 K) and for the metals with high-temperature martensite transitions which 
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are usually accompanied by essential increase in the AE in the lattice dynamics [fL0|j . In 
the present work the AE features (in the comparison with the electron contributions) in 
the thermodynamical properties of refractory fee metals Ir and Rh and in the bec and 
fee phases of Ca and Sr near the points of structural transformation are studied. The 



calculations were made in terms of microscopic models used in |10|, [TTJ] which describe a 
wide range of the lattice properties of these metals. 
The initial Hamiltonian of phonon subsystem is 

A V 7 n=3 

Ai ...A n 

where M is the mass of atoms; uj\ is the phonon frequency, A = q£ where q is the wave 
vector; £ is the branch number, Q\, P\ are the normal coordinates in the harmonic ap- 
proximation and the corresponding pulses, are amplitudes of intermode interactions. 
In the leading order of the parameter 77 = T/E at where E at is the energy of the order of 
binding energy, the basic contribution of AE to free energy F is determined by the second 
order by and the first order by H^> [p], [J 

F an = F® + F (4) (3) 

where at T Q D (Q D is the Debye temperature), 

T 2 V ^ (3) (W)| 2 () 

F (4) T 2 ^ $ (4) (A, X,fi,fi) _ 

X/j, A ^ 

as usually, summing over phonon quasimomenta is fulfilled taking into account the conser- 
vation law for A = (£, k), fi = (77, q), = (C, k + q). The explicit microscopic expressions 
for $® and $< 4 > are given m y|. 

The heat capacity of the crystal at T > S^, without the standard harmonic contri- 
bution C v arm = 3R (where R is the gaseous constant), has the form: 

AC = C P - C v arrn = (C P - C v ) + C v n + C e v + C^ (6) 

where Cy n are the anharmonic contributions to the heat capacity 

00 

df(Ey 



C V (T) = I J dE 
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(E - ^N(E) (7) 



is the electron heat capacity (N(E) is the electron density of states, is the chemical 
potential, f(E) is the Fermi distribution function), C v is the contribution of defects; 
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the calculation formula for Cp — Cy in the quasiharmonic approximation is presented in 
@. The density of states N(E) was calculated by the FP-LMTO method in the local 
approximation for the density functional (LDA) ||12j| . To estimate the contribution of 
defects to the heat capacity a common approximation of independent monovacancies was 
used with 



ct 



(T) 



cxp 



(8) 



TT|, show that near the 



where E v , S v are the energy and entropy of monovacancy formation, respectively. 

The calculations of AE in the Ir lattice dynamics, made in 
melting temperature T = T m the temperature dependences of the phonon frequency and 
phonon damping drastically change over the Brillouin zone reaching the values of the order 
of 15-20%. This value is high enough so that the question on the role of higher orders in 
the anharmonic perturbation theory could be set. The theory of self-consistent phonons, 
often used for this purpose (see, e.g., [|U|]) is not adequate in this case as it only accounts 
for a part of contributions to ([!]) and neglects the if( 2n+1 ) contributions. At the same 



time, the calculations ]Tl[ show that in Ir the contributions of three-phonon processes to 
the temperature shifts of the phonon frequencies always dominate over the contribution of 
four-phonon processes. In the recent work ||] another approximation is proposed which 
accounts both the three-phonon and four-phonon processes, and, as stated, reasonably 
describes the thermodynamics of crystals with the Lennard- Jones type interaction. The 
applicability of this approximation in the case considered seems somewhat doubtful as 
it is based on the replacement of anharmonic frequency shift by a certain value average 



over the Brillouin zone while the calculations in |]TT| show that the q-dependence of AE is 
essential. Moreover, it is known that the approximations adequate for the description of 
thermodynamics of systems with "hard" potentials of Lennard- Jones type are not usually 
applicable for metals as they have relatively "soft" potentials ||. 

To calculate the effects of higher orders by the Hamiltonian the approximation 
similar to that of Gell-Mann and Brueckner ("bubble approximation") in the theory 
of electron liquid |TJJ was used. It can be shown that the accurate expression for the 
contribution of three-phonon processes to the free energy at T ^> Bd has the form: 



^ (3) = -|/^E^(- = o)s^ = o)] 



(9) 



where G\, H\ are the Green function and self-energy part for phonon (see |l 
by the Dyson equation 
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The calculations [|11[| show that the AE in fee metals have pronounsed peaks as functions 
of q at points X and L. The largest contribution to the free energy comes from the 
"bubble" diagrams; this situation is similar to the electron gas of high density with the 
only difference that in the latter case the "dangerous" region is q ~ |14]]. Taking this 
into account, in ([5|)-(|n]) we may restrict ourselves to the expression of the lowest (second) 
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order of the perturbation theory for Then we have 
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Unlike the bcc metals the fee metals have no soft modes in the phonon spectrum, which 
makes the four-phonon processes relatively small comparing with the three-phonon ones 
|10| , |TT|1 . Therefore for it is sufficient to use the simplest corresponding to the Hartry- 
Fock approximation (the simplest procedure of disengagement in averaging Hamiltonian 
H^). Then for T » 0^ we have 
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where 
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is the corresponding correction to the phonon frequency. 
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Fig. 1. Anharmonic contributions to the free energy of Ir and Rh (in K) considered by the perturbation 
theory (fj), (^) ("old") and with taking into account the higher orders (|TT1), ( ^2| ) ("new"). 



The computational results are shown in Figs. 1-5. One can see that the AE higher order 
contributions to the free energy of Ir and Rh are not large (do not exceed 20%). In Ir and 
Rh the anharmonic contribution to heat capacity is small comparing with the electron 
contribution while in the sp-metals (as the hypothetical fee phase K shown in Fig. 3) 
these contributions are comparable. Among the metals studied only for Ir the detailed 
experimental information is available, the theoretical data being in good agreement with 
these results (Fig. 4). In the estimation of Cy according to (|8|) the the first principles 
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Fig. 2. Electron contribution to the heat capacity (in R units) and temperature dependence of chemical 
potential (to 1CP 4 Ry) in Ir; the dependencies for Rh are similar. 
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Fig. 3. Comparison of the electron and anharmonic contributions to the heat capacity (in R units) for 
the fee phases of K and Ir; t — T/T m . 
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Fig. 5. Theoretical value of heat-capacity in Ca. The notation introduced C v = C^ arm + C°" + Cf, 

s~lharm i /^tan i riel i ( /~i f~i \ r~if riharra i r^an i r^free el r^f r^harrn i /~<an i r^free el i 1 s~l f i \ 

~ +0„ +U V +(Lyp — Ly v ), Ly^, — L- v +<~' v , — L- v +<^ v +[Ly p — 0„J, 

(jfree el _ fog electron heat capacity in the approximation of free electrons. 
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calculation result from E v = 1.88 eV [I7| was used for E v , and S v = 2.5 for S v 0. It 
follows from Fig. 5 that in the difference of heat-capacity for bcc and fee phases of Ca (and 
similarly Sr) a strong compensation for different contributions - electron, anharmonic ones 
- and the contribution of Cp — Cv takes place. Thus, for the discussion of thermodynamics 
of high temperature structure transformations in metals the microscopic calculations of 
corresponding values including the anharmonic effects are needed. 
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